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Abstract 


In this paper, we use a novel technique to solve the nonlinear fractional 
Volterra-Fredholm integro-differential equations (FVFIDEs). To this end, 
the Legendre wavelets are used in conjunction with the quadrature rule for 
converting the problem into a linear or nonlinear system of algebraic equa- 
tions, which can be easily solved by applying mathematical programming 
techniques. Only a small number of Legendre wavelets are needed to obtain 
a satisfactory result. Better accuracies are also achieved within the method 
by increasing the number of polynomials. Furthermore, the existence and 
uniqueness of the solution are proved by preparing some theorems and 
lemmas. Also, error estimation and convergence analyses are given for the 
considered problem and the method. Moreover, some examples are pre- 
sented and their results are compared to the results of Chebyshev wavelet, 
Nystrom, and Newton—Kantorovitch methods to show the capability and 
validity of this scheme. 
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1 Introduction 


Fractional calculus has been applied in physics in recent years, although it has 
a long history in mathematics. Most of the physical phenomenon can be mod- 
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eled by using fractional calculus [34, 10]. Applications of fractional differen- 
tial equations and fractional integral equations create a wide area of research 
for many researchers [3, 4, 12, 11]. It has been applied to model the non- 
linear oscillation of earthquakes, fluid-dynamic traffic, frequency-dependent 
damping behavior of many viscoelastic materials, continuum and statistical 
mechanics, solid mechanics, economics, signal processing, and control theory. 
To better analyze these systems, it is required to obtain the solution of these 
equations. However, mostly, finding the analytical solution of these equa- 
tions cannot be possible, so considering more accurate numerical solutions 
can be helpful. In the literature, there are various techniques for solving 
fractional ordinary differential equations (FODEs), fractional partial differ- 
ential equations (FPDEs), fractional integro-differential equations (FIDEs), 
and dynamic systems with fractional derivatives, such as analytical and semi- 
analytical methods (homotopy analysis method, Adomian’s decomposition 
method, etc.) [1, 13, 14, 15, 17, 39, 35] and numerical methods (finite differ- 
ence schemes, collocation method, Tau method, etc.) [6, 36, 16, 19, 30, 5, 22]. 


We can find some other famous methods for the numerical solutions of 
these kinds of equations. For example, spline collocation [29], analytical 
Lie group approach [28], fractional differential transform [25], least-squares 
[20], rationalized Haar function [27], exp-function method [8], and many oth- 
ers. In [32], a novel Legendre wavelet Petrov-Galerkin method was pre- 
sented for fractional Volterra integro-differential equations. The Chebyshev 
wavelet method [23] has been used to nonlinear fractional Volterra-Fredholm 
integro-differential equations (FVFIDEs) with mixed boundary conditions. 
Approximate solutions of Volterra-Fredholm integro-differential equations of 
fractional order based on the Sinc-collocation method have been discussed 
in [2]. In [24], the Nystrém and Newton—Kantorovitch methods were de- 
scribed for solving FVFIDEs with mixed boundary conditions. Systems 
of integro-differential equations [37] have been solved by using the spectral 
second kind Chebyshev wavelets scheme. Also, nonlinear Volterra integro- 
differential equations of fractional order have been solved numerically using 
the Euler wavelet method [38]. The reliable modification of the Laplace Ado- 
mian decomposition method [9] has been applied to solve nonlinear Volterra- 
Fredholm integral equations. Haar wavelet bases have been employed for solv- 
ing the integro-differential equation [7]. In [18], hybrid Bernstein block-pulse 
functions have been used for solving systems of fractional integro-differential 
equations. 


Orthogonal polynomials and functions apply to different problems be- 
cause of their appropriate attributes. These functions and methods such as 
collocation, Galerkin, and Tau, are applied to reduce the solutions of different 
problems to the solutions of a system of algebraic equations. In this work, 
the Legendre wavelets are implemented to obtain an approximate solution 
for nonlinear FVFIDEs. We consider the fractional FVFIDEs in a Banach 
space as follows: 
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(9 Doy) (e) + Vy(z) + Fy(x) = g(a), m-1<agm, meN, (1) 


where § D? is in the Caputo sense fractional derivative, g(x) is a given con- 
tinuous function, x € J = [a,b], y(a) is unknown function, and 


‘eS b 
Volz) = fhe, 9)Na(ule)) dé, Pye) = bole, )Nalu(e)) a8, 

in which k; : Jx J > R,i = 1, 2, are continuous functions and N;: JxR- R, 

i = 1,2 are Lipschitz nonlinear continuous functions. Equation (1) is subject 

to the following mixed boundary conditions: 


m 


S- sy? @) + my) (b)| =, 1=1,2,...,m. (2) 


=1 


The main characteristic of this technique is that it reduces these prob- 
lems to a system of algebraic equations. This approach is based on converting 
the FVFIDEs into mixed Volterra-Fredholm integral equations through inte- 
gration, approximating various terms involved in the equation by truncated 
Legendre wavelet series and using the operational matrices to eliminate the 
integral, derivation, and product operations. 

The rest of the paper is organized as follows: In section 2, some es- 
sential mathematical preliminaries and definitions of fractional calculus are 
introduced. In section 3, the properties of Legendre wavelets and function 
approximations are discussed. In the next section, the Gauss quadrature 
Legendre wavelet method (GQLWM) is constructed for solving FVFIDEs. 
In section 5, we study the convergence and error analysis of our algorithm. 
The existence and uniqueness of the solution are investigated in section 6. In 
section 7, the proposed method is applied on two examples to demonstrate 
the efficiency and accuracy of the present method. At last, a brief conclusion 
is given in section 8. 


2 Basic definitions and notations of the fractional 
calculus 


In this section, some definitions and properties of the fractional calculus, 
which will be used, are presented. For more details, see [21, 26, 31, 33]. 


Definition 1. The Gamma function is intrinsically tied with fractional cal- 
culus. The definition of the gamma function is given by 


T(a) = I e §€% dé, Re(a) > 0. 
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Definition 2. A real function g(x), x > 0, is said to be in the space Co, 
6 € Rif there exists a real number p > @ such that g(x) = x?g;(x), where 
gi(x) € C (0,00) and it is said to be in the space Cf, k € NU{0} if and only 
if g”) € Co. 

Definition 3. The Riemann-Liouville fractional integral operator of order 
a > 0 of a function g € Co, 0 > —1, is defined as 


ee) = ao | *(@ — 1 9(é)dE. 


Some properties of the Riemann-Liouville fractional integral are as fol- 
lows: 


1. of?g(x) = g(z), 
2. (ol o129) (x) = (ofB ol) (x) = (oft 9) (2), 


3. gI@ (a — p)* = ee py, SO Ae Me 


Similar to the integer-order integration, the Riemann-Liouville fractional in- 
tegral is a linear operator, which means that 


oly (>. afte) = be -) oft fil), 


i=1 


where {c;}’~, are constants. 
Definition 4. The fractional derivative of g(x) in the Caputo sense is defined 
as 

0 Deg(2) = ofp “9 (a), 
forn-—l<a<n,nEN, x >0, andge C”,. 


The relation between the Riemann-Liouville operator and Caputo opera- 
tor is given by the following lemma. 


Lemma 1. Ifn-—1<a<n,né€N, then § D%9I”~“g(x) = g(x) and 


n-1 k 
m—Qa Qa x 
olf 6 De g(x) = g(z)— Yi gM OT, «>. 
k=0 


3 Legendre wavelets and functions approximation 


3.1 Properties of Legendre wavelets 


Wavelets are a family of functions constructed from the dilation and trans- 
lation of a single function called the mother wavelet. When the dilation 


IJNAO, Vol. 12, No. 1, (2022), pp 229-249 


Legendre wavelet method combined with the Gauss quadrature rule for ... 233 


parameter a and the translation parameter b vary continuously, we have the 
following family of continuous wavelets as 


b 
Vap(z) = |al72U (: = -) , a, DER, a¥0. 


If we restrict the parameters a and b to discrete values as a = ag kb= 
nboag eg 0, bb > 0 and n and k are positive integers, then the following 
family of discrete wavelets will be obtained 


brn (@) = |ao|? (aka — nbo), 


where W%.n(x) form a wavelet basis for L?(R). In particular, with ag = 2 and 
bo = 1, Yen(x) form an orthonormal basis. 

Legendre wavelets Wn (a) = (k, i,m, az) have four arguments; 2 = 2n — 1 
for n = 1,2,...,2*-!, k € Z+, mas the order of Legendre polynomials and 
x, which is the normalized time. They are defined on the interval [0,1) by 


2°= (2m + 1)? L* @), *2<e¢<22 
Wnm(£) = W(k, n,m, x) = m ; S 


0, otherwise, 


where m = 0,1,..., M—1, n = 1,2,...,2*~1, the dilation and translation 
parameters are a = 2~* and b = f2-*, respectively. Moreover, L* (x) = 
Lm(2"x — ft) are shifted Legendre polynomials on interval [0,1), which are 
orthogonal with respect to the weight function w(z) = 1. Also, L,,’s can be 
determined by the following recurrence formulas: 


3.2 Function approximation by Legendre wavelets 


A function y(x), which is square integrable in [0,1), may be expressed in 
terms of the Legendre wavelet as 


a=) YS triPaanle), (4) 


n€N meENU{0} 


where Ynjm = (y(£),Wn,m(x)) and <-,- > denotes the inner product. If the 
infinite series in (4) is truncated, then one obtains 
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2k-1 M1 


y(t) ~ S~> S~ Yn Pnm(#) = YF U(x), (5) 


n=1 m=0 


where the coefficient vector Y and the Legendre wavelet function vector V(x) 
are given by 


T 
Y = [y1,0, Y1,15 +++) Y1,M—1) Y2,05 Y2,15 +++) Y2,M—1y+ +63 Yok-19)Yok—-1,45+++sYok-1, il > 
(6) 


T 


W(x) = [v1,0(x), y11(2), trey W1,M-1(2), ae »Wok—1 (2), Wor—11(2), tee »Wok-1 yg—1(2)] : 
(7) 


Similarly, we can approximate the functions k;(xz,€) € L?((0,1] x [0,1]) as 
k;(a,€) = U(x)" K,W(E), (8) 


where K;, i = 1, 2 are 2*-!'M x 2*-!M matrices. 


4 GQLWM for solving FVFIDEs 


Let us consider the nonlinear FVFIDE with mixed boundary conditions given 
by (1)-(2). To approximate the functions y(x) and g(x) by y(x) = Y7U(z) 
and g(x) = G7 W(x), assume that 


Ni(y()) = ul), No(w(€)) = v(§), (9) 
where u(£) and v(€) are as follows: 
u(£) =U WE), o(€) =V7 VEE), (10) 


in which U7 and V" are defined similarly as in (6). Applying the operator 
of@ on both sides of (1), results in 


= es ak 1 ie a—1 
10) — DOF TG [ @-1) oar 
1 . a-l r 
Pa (x —T) i: ky (7, €)u(§)drdé 
1 . a-l : 
Sal (= 5) i ko(r, €)u(€)drdé. 


Replacing the exact solution y(a) by Y7W(z) and using their approximations 
by (5), (8)-(10), one gets 
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n-1 


YT W(x) ye (0+) a 
=Tay f, Mr 
— wy f @- 7 [ vo" KweuTw@arag 
7 Fa) i ae [ U(r)" KoW(E)VTW(E)drdé. (11) 


Collocating (11) in M2*~! Legendre wavelet-Gauss collocation points 2;, 
j =1, 2,..., M2*-!, on interval [0,1], we arrive at the following system: 


n—-1 ak 
YT W(2;)- >> se aE 
k=0 


1 a) ne 
=a I, (aj — 7)°-1G?W(r)dr 
7 aot J (xj — eae U(r)? Ky W(E)UT W(€)drdé 
= ns i (aj —7)*"" [ U(r)? Ky W(€)VEW(E)drd€. (12) 


In a similar way, the mixed boundary conditions (2) are approximated as 
follows: 


‘ [u¥? YI (0) = ng¥TWO-Y(1)) =, i=1,2,...,.m. (13) 


j=1 


Now, we use the quadrature rule to approximate the integral involved in 
this equation, which has zero error integration for polynomial integrands of 
degree less than or equal to 1 + M2* with M2*-! Legendre-Gauss nodes. 
To do this, M2*~! intervals [0,;] are transferred to a fixed interval [—1, 1] 
through the transformations tT = (x;/2 )(s+ 1). Applying the Gaussian 
quadrature, system (12) is converted to 


n-1 gk 
¥?U(2;)- >> ery 


k=0 
1 ee) ae 
=r /, (xj — 7) GE W(r)dr 
1 M2k-1 a sana tie 
Ka? 2 wi(F- 9) [wort evequ™ eas 
1 ne 2s a-1 ft 
~ Fey 2 Py wi(Fa-s) [ wert wovrweus, (14) 
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where o = *4(s;+1)) and w for! = 1,2,...,M@2*~" are the weight functions 
on [—1,1]. Taking into account (9)-(10) and using collocation points 2,, 
j =1,2,...,M2*-1, one obtains 


Ny(¥TW(2;)) =UTW(a;), No(YTU(e;)) =VTW(e;). (15) 


Combining (13)—(15), the main problem is reduced to a system of nonlinear 
algebraic equations. By solving this system, the unknown values Y7 will be 
obtained, and by relation (5), the approximate solution of system (1) will be 
determined. 


5 Convergence and error analysis 


In this section, we will obtain an estimation for the error bound of our numeri- 
cal method. Also, its convergence analysis will be discussed. For this purpose, 
assume that L*, (2) = Lm(22—1) and that A(x) = [L*(x), L*(x),..., Ly (2)|". 
So, the function y(x) € L?[0,1] can be expressed in terms of the Legendre 
polynomials basis A(z) as 


N 
yo) = So nbs (x) = LTA(a), 
n=0 


where L = [I,l,...,ly]°. 


Lemma 2. Suppose that y € C™*?(0, 1] and that S = span { Lj, Lj,..., Ly} 
is a vector space. If L7 A(z) is the best approximation of y(x) out of S, then 
the error bound of presented method is as follows: 


2y2™43 Mint 
m+ Ih/2m +3’ 


in which M441 = max {|p "(@)| : O<a2<1} and y = max{1 —€,€}. 


BR 


lu 2TAI, < 


Proof. The Taylor polynomials, implies that 


(x = ey ” 


I y (6) 4. ¢ CED ym 


mi 


Ty(x) = yo(€) + (w — €)y"(€) + (6), 


where € € (0,1). Therefore, there exists \ € (0,1) such that 
(x = ie m 
Iu(x) ~ Ty(x)| < |{F—P— yr (ay). 
Since L7 A(x) is the best approximation of y(2), one obtains 
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1 
lly — 7A]? < Iv Tul? = i. ly() — Ty(2) 2dr 


1 _ m+ 
< | 
0 
2 


Degree ss 
~ (2m + 3)[(m +1)!’ 


and the proof is completed. 


Theorem 1. Assume that y € C™**[0,1] and that its approximation by 
the Legendre wavelets is 9(a) = Y7 W(x). Then, its mean error bound is as 
follows: 
oer) (8 —§) 4/9 Mi 
Yl, < 
lu tla S Ga Dam +3 


Proof. Dividing [0,1] to 2*-! subintervals Ix.» = [(n — 1)/2*-1! ,n/2*-? J, 
n = 1,2,...,2*-1 with the restriction that 7 is a polynomial of degree less 
than m+ 1 and also using Lemma 2, we get 


ie as fm ly(a) — G(a) Pde 


=> [P Wye)- g(a) P ae 


(a= WGrets) 9(1—k)(2m-+3)+1 gk-1 


~ [m+ )IPFem+3) & 


J2M,,2 
<> (m+ 1)h/2m +3 


iments 4 


[(m + 1)! (2m +3) ’ 


where M,, = max {|y(™+) (2) 


, cel, i which completes the proof. 


Lemma 3. Suppose that for y(x) € [0,1), there exists 8, € R such that 
ly(x)| < By. Then, the sum of Legendre coefficients of y(x) defined in (5) is 


absolutely convergent if |Ynm| << V2!-*By. 


Proof. The function y(z) € L[0, 1] can be expressed as the Legendre wavelet 
basis as defined in (5). For m > 0, we have 
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nl =1<t¥om >1=|f uf ‘i. asl S < [iw ‘ifieastlee 


< By [ i @leak | its sae) 


nk 


= 5,2°F* am +14 f |Lin(2*a — 2n + 1)| de. 


Ink 


Setting s = 2*2 — 2n +1, one obtains 


= (2m +1)2 / |Lm(s)|ds. 


-1 


[Yn,m| < By2 


The Holder’s inequality implies that 


1 
|Yn,m| < By2~ “Fam +1)! = vr * By. 
2k~1 M-1 
This means that the series S> 3° Yn.m is absolutely convergent as k > 
n=1 m= 


ow. 


k-1 
Theorem 2. The series solution §(#) = S77 _, y ey mWn,m(«) defined 
by (5) is convergent with respect to L*- “norm on (0, 1] if the sum of absolute 


M-1 
values of the Legendre coefficients see ‘i i cue 6 [Ynm| for continuous function 
y(x) forms a convergent series. 


Proof. Let L*(IR) be a Hilbert space and let Wp defined in (3) form an 
orthogonal basis with W(x) = Aj (x) and 6; =< g(x), A;(x) >. Define a 
sequence of partial sums {,S,,} as 


x) = > 7 5;A;(2) 
j=0 
Now, for every € > 0, there exists N. > 0 such that for every n > m> N, 
[XS sane 
l=m-+1 
= do 1a. 


l=m+1 


Sn(2) — Sim(x) [2 dx< S|? [ \Ai(2) Pde 


l=m+1 


Using Lemma 3, 577", |61|? is absolutely convergent. Hence, according to the 
Cauchy criterion, we have 
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n 


S- \o1|7 < : 


l=m+1 


So, ||Sn() — Sm(x)||, < Ve? =e. Thus, the sequence of partial sum of the 
series is a Cauchy sequence. Therefore, it convergent with respect to L?-norm 
and this completes the proof. 


6 Existence and uniqueness 


In order to study the uniqueness of the solution, we consider the FVFIDEs 
(1) in an operator form as 


(§ Dey) (a) = g(a) — KiNiy — K2Noy, (16) 


where 


z b 
KiNiy =f ha(e.8)Ni(y@)) as, KaNay = f ba(2,6)Na(ul6)) 4. 
Applying the operator oJ on both sides of (16), one obtains 
y(x) = f(x) + off (g(x) — KiNiy — K2Nay), 


where f(x) = y) (0+)a*/k!. This equation can be written as Ty = y, 
where 7 is as follows: 


Ty(x) = f(x) + off [g(x) — Ki Niy — K2Noy]. 


Let (C'[0, 1], |oo||) be a Banach space of all continuous functions with the 
norm ||h||oo = max,;|h(x)|. Also, the operators N; and N2 satisfy the Lips- 
chitz condition on [0,1] with Lipschitz constants J; and Jz as 


> |Noy(x) — Nay(x)| < Jo |9(@) — y(a)]. 


INig(x) — Miy(x)| S$ Ai |9(@) — y(@) 


In the following theorem with these assumptions, we show that the FV- 
FIDEs (1) have a unique solution. 


Theorem 3. Suppose that the nonlinear operators N; and No satisfy the 
following relation 


Fill Nilloo + Fall Nello < F(a + 1). 
Then the FVFIDE (1) has a unique solution. 


Proof. Let T : C [0,1] + C [0,1] be defined as 
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Tule) = $2) + seas fe" 6) ~ KiNi(@) — Kaylee 

Suppose that gy, y € C'[0,1]. Then for every positive x, we have 

Ty(x) — Ty(z) 

= maf (x — €£)°"* [Ki Niy(€) — KiN19(€) + KoNoy(€) — K2Nog(é)] dé. 
Pa) Jo 

Therefore, one obtain 

IT g(x) — Ty(2)| 


1 . a—1 de - 
Ta) ‘ le — €°* [|| [Niy(€) — Nag(©)] + [Kal [Nay(€) — Nog(6)l] dé 


< Fay A la — €[-* [|B] Ja |G(€) — (| + [Kol Je [() — wlll a 


1 ° a-1 2: 
< rai. Je = €°* [ Kaloo + Kalo Jallld(€) — yO)lloo 4 
= NCES : iy Ml Kallac Ft + Ilo 42) IE) — Yc 
< Tea (| illoodt + Belloc J2) IGE) — y(8)lloo- 


Hence, 


ITG(x) — Ty(@)lloo S HH) — y(lloo 


where L = [1/[(a +1) ] (||Ai||,,41 + ||Kal|,,J2). Since L < 1, by the con- 
traction mapping theorem, this problem has a unique solution in C (0, 1]. 


7 Numerical examples 


In this section, to demonstrate the capability and accuracy of our method, we 
apply it to some examples and compare the quality of the computed solutions 
with some other well-known methods. In order to demonstrate the error, we 
will calculate the maximum absolute error (MAE) by the following formula: 


MAE = mex {lul@) =gle)lt 


Also for comparison, in Example 2, the root-mean-square error (RMSE) is 
calculated. This error is as follows: 


RMSE = i i [y(x) — g(a)|Pde. 
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Example 1. Consider the following nonlinear FVFIDEs: 


(Spey) @— [AEE ac [c+ 20 M1608 = a0, 


subject to the boundary conditions y(0) = 0 and y(1) = 2 in which 


The exact solution of this problem is y(a) = x? +a. Numerical results for 
Example 1 are indicated in Table 1, which show that the results of our method 
are better than the results obtained in [24, 23] by the Chebyshev wavelet 
method (CWM), Nystrom and Newton—Kantorovitch methods, respectively. 
Table 2 analyzes the exact solution y and the approximate solutions of this 
algorithm by CWM and GQLWM. One can see that the present method has 
excellent accuracy with respect to CWM. Table 3 shows the maximum ab- 
solute error between the exact and approximate solutions for various choices 
of M and k. The outcomes reveal that the results of GQLWM by using 
only a small number of bases are very promising and superior to CWM and 
Nystrom and Newton—Kantorovitch methods. Furthermore, we compare the 
absolute error functions for different values of k and M in Figure 1. These 
results confirm that with increasing the amounts of k and M, the error will 


be decreased. 
YY 
10-2 _ a wrest ety, oe 


oe eg are es 


" 
oe ww, oe ee Oe 

wif : : . 7 : ; He 
10°" . 


\ Bhan i ' 
rr t 1 Y aa aul ‘an my mv 


. ‘i ry eee s 2 hry parent 


0 02 04 0.6 08 i 


— —* 
— 


— M=4= = 1 M=6 = = M=8 =" M=10 — M=4= = 1 M=6 = = M=8 =" 'M=10 


(a) k=1 (b) k=2 


Figure 1: Error history of the presented method for various values for M and k of 
Example 1. 
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Table 1: Comparison between our method, CWM [23] and method in [24] for Example 
1p 


Ma 0.1 0.3 0.5 0.7 0.9 
k=3,n=160 | 781E—-06 155F-05 165EF-05 127E-05 5.07E—-06 
Error [24] k=3, #7 = 320 192E—06 3.84F-06 410F—-06 3.15E-—06 1.25E—06 
k=3,n=640 | 4652-07 931E-07 9.92E-—07 765E-—07 3.04E—07 
Error [23] fork =1,M=19 | 347E—-12 608E-12 621H—-12 468H-12 1.84k—12 
k=2, M=4 2.93-—06 1.01f-—05 1.99F-06 1.08h-05 4.27 — 06 
k=3, M=5 | 817E-09 5112-10 505E-11 6.78E-10 9.48E -—09 
GQLWM k=4, M=6 169F —23 4.22F-—20 991F—-—21 7.06F-—20 4.45F — 23 
k=1,M=10 101F-12 467£-13 100H-12 4.71F-13 1.02E-—12 
-=1,if7 = 19 L16L-28 211F-30 886£-—29 3.35E-—30 1.30E — 28 
Table 2: The numerical solution of Example 1. 
x | y(exact) y(CWM) y(GQLWM) 
0.0 0.00 0.00000000000000000001 0.00000000000000000000 
0.1 0.11 0.11000000000347509219 0.10999999999999999999 
0.2 0.24 0.24000000000518991700 0.23999999999999999999 
0.3 0.39 0.39000000000608906284  0.38999999999999999999 
0.4 0.56 0.56000000000639391977 0.55999999999999999999 
0.5 0.75 0.75000000000621724965 0.74999999999999999999 
0.6 0.96 0.96000000000562994250 0.95999999999999999999 
0.7 1.19 1.19000000000468149700  1.18999999999999999999 
0.8 1.44 1.44000000000340880150  1.48999999999999999999 
0.9 1.77 1.71000000000184060200 = 1.70999999999999999999 
1 2.00 1.99999999999999999990 = 2.00000000000000000000 


Table 3: Maximum absolute error between the exact and approximate solutions for 


various choices of M and k for Example 1. 


M | 4 


7 10 


13 16 19 


k=1 
k=2 
k=3 


441E—-04 4062-08 1.03F-12 402h—-16 3.01E — 23 
768E—-05 941EF-10 3052-15 3968-21 2.55 — 27 
5.46£-—06 8.23F—-12 3398-18 549F-25 4.42F — 32 


1.57E — 28 
1.61LE — 34 
2.02E — 39 


MAE 


Example 2. Consider the following nonlinear fractional integro-differential 


equation: 
1 
(¢ Dey) (2) = g(x) + | xe 
0 


where 0 < a < 1 and g(z) is as follows: 


[y(é))’dé, y(0) = 0, 


1 


Case 1. [23] g(x) 4 


y(x) =a. 


= sry tet 3 


3 Z and exact the solution is 


x with a = 
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Case 2. joes Sr (2) V2 gt ;@ with a = } and the exact solution is 
y(vz) = 2°. 


Case 3. [40] g(x) = 1 — $2 with a = 1 and the exact solution is y(x) = a. 


We apply our method to this example with various values of M and k, 
and the results are presented in Tables 4 and 5. The maximum absolute error 
of y(a) by GQLWM for k = 1 with different values of M has been compared 
to the error of CWM [23] for cases 1 and 2. These results demonstrate the 
validity and capability of GQLWM with respect to CWM. The graphs of 
absolute error functions for different values of M and k with a = 1/4, a =1, 
and a = 3/4 are given in Figures 2, 3, and 4, respectively. These figures reveal 
that the error will be decreased when M is increased. Table 5 compares the 
RMSE of GQLWM with the results in [40] for different values of k and M 
when a = 1. This table reveals that, for a certain value of k, as M increases, 
the accuracy increases. Also, for a certain value of M, as k increases, the 
accuracy increases as well. Therefore, GQLWM for solving this problem is 
very effective and more accurate with respect to CWM and the method in 
[40]. 


Table 5: RMSE for some k and M of GQLWM and the second kind Chebyshev wavelet 


method [40] for Example 2 (case 3). 
Method in [40] GQLWM 
M 2 2 6 10 14 
k=3 2.970E — 07 2.348F —09 4.2368-—11 1.003F—12 1.320F — 13 
k=4 1.861E — 08 7.391F—-10 9.86838 —-13 4.002F—15 7.5691E —17 
k=5 1.164E — 09 8.541F—12 8.786£—-16 9.251F—-18 2.015£ — 20 


(a) k=1 (b) k=2 


Figure 2: Error history of the presented method for various values for M and k of 
1 


Example 2 with a= 3. 
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Table 4: Comparison between our method and CWM [23] for Example 2 (cases 1 and 2) when k = 1. 


a a= zi a= a 
Error[23] GQLWM Error[23] GQLWM 

x M=19 M=13 M = 16 M=19 M=19 M=13 M=16 M=19 
0.0 | 1.68 — 20 0 0 0 7.92E — 21 0 0 0 
0.1 | 906F-—13 175kH-13 185F-16 1.27F-20 | 6.57E—-18 653EF-16 400F—-19 2.51F —24 
0.2 |399F-12 4.05F-12 7538-17 3.92F—23 | 362F—-18 589F-15 1683H-19 2.91E—25 
0.3 | LOBE-11 3.76F-12 147F-16 149F —20 | 265F—-18 544F-15 3.80F-19 6.18E — 24 
0.4) 2.10hF-11 2.21F—-12 1.538£-16 1.068 —20 | 2.11kF-18 341£-15 4.19F-19 4.08F — 24 
0.5) 3.73F-11 835F-14 148F-16 1838-22] 185F-18 407E—-16 425F-—-19 2.29F — 25 
0.6) 605F—-11 238F-12 142F-16 104f-—20]165F-18 261F—-15 428F—-19 4.33B — 24 
0.7 | 916-11 3.92F—-12 1.23F-16 1518-20] 149F-18 463F—-15 3.94EF—-19 5.92KF —24 
0.8} 131F-10 423F-12 148f-16 3.71F —22 | 1.40£-—18 7 L8lE-—19 1.23E — 25 
0.9} 182F—-10 3.11F-13 167E-17 1.32F-—20 | 107E—-18 498F-15 510EF-19 4.92E —24 
1 | 2.452 — 10 0 0 0 5.42E — 17 0 0 0 
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(a)k=1 (b) k=2 


Figure 3: Error history of the presented method for various values for M and k of 
Example 2 with a = 3. 
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Figure 4: Error history of the presented method for various values for M and k of 
Example 2 with a = 1. 


8 Conclusion 


In this paper, the GQLWM, which is an efficient technique for the numerical 
solution of nonlinear FVFIDEs, has been proposed. Using this method, the 
system of nonlinear FVFIDEs was reduced to a system of algebraic equa- 
tions. The numerical solution of the resulted system was approximated by 
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the Gauss quadrature formula with respect to the Legendre weight func- 
tion. By solving this nonlinear system, the numerical solution was obtained. 
Moreover, the convergence, error analysis, existence, and uniqueness of the 
proposed method were discussed. It was evident that, for a certain value of 
k, as M increases, the accuracy of the GQLWM is increased. Also, for a 
certain value of M, as k increases, the accuracy is increased, as well. The 
method was applied to two examples, and the obtained results were com- 
pared to some other well-known methods. This comparison showed that the 
GQLWM is a suitable and powerful technique for solving the nonlinear FV- 
FIDEs. In the end, we note that the method can be easily extended and 
applied to multi-dimensional integral equations or systems of FVFIDEs eas- 
ily with some modifications. We also believe that it shall not be difficult to 
extend this approach to nonlinear equations of general form, which will be 
the subject of future researches. 
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